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Abstract. The effect of shear flow on the phase-ordering dynamics of a binary mixture with field-dependent 
mobility is investigated. The problem is addressed in the context of the time-dependent Ginzburg-Landau 
equation with an external velocity term, studied in self-consistent approximation. Assuming a scaling 
ansatz for the structure factor, the asymptotic behavior of the observables in the scaling regime can be 
analytically calculated. All the observables show log-time periodic oscillations which we interpret as due to 
a cyclical mechanism of stretching and break-up of domains. These oscillations are dumped as consequence 
of the vanishing of the mobility in the bulk phase. 

PACS. 47.20Hw Morphological instability; phase changes - 05.70Ln Nonequlibrium and irreversible ther- 
modynamics - 83. 50 Ax Steady shear flows 

1 Introduction length scale, which grows as a power law L(t) ~ t z . The 

spatial patterns of the domains at two different times are 

A well-studied problem in non-equilibrium statistical me- related b y a S lobal cnan S e of this len g th scale - A si g na " 

chanics is the growth of domains in quenching processes ture of this dynamical scaling is the fact that the structure 

0. Namely, when a system is suddenly quenched from a f actor C(k, t) can be cast in the form 
disordered initial state into a thermodynamic region where 

C(k,t) = L d (t)C (kL(t)) (1) 

different phases coexist, macroscopic domains can be ob- 
served, usually characterized by a single time-dependent where C Q is the scaling function. 
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In the case of a binary mixture, the evolution of the 
system is described by a scalar order parameter (j) rep- 
resenting the difference of concentrations of the two liq- 
uids in the mixture. When hydrodynamic effects are not 
considered, the quenching process can be described by an 
equation of the form d t (f> = — Vj, where j = —rV(SH/S(p), 
H is a free-energy functional describing the ordered phases 
and r is a transport coefficient called mobility, usually 
taken constant. It has been argued [0 that the dynamics 
is more accurately mimicked considering a field-dependent 
mobility of the form 

m = (l-a(T)^) (2) 

where a(T) — > 1 for temperature T — > 0, a(T) —* for 
T — * T c . In this way it is possible to take into account 
the different nature of the two mechanisms which operate 
during the phase separation: surface diffusion and bulk dif- 
fusion. The first has a growth exponent z = 1/4 || and is 
due to the diffusion of molecules of the two species along 
the interfaces. The second, which is also called Lifshitz- 
Slyozov mechanism, is due to the diffusion of molecules of 
one species from more curved interfaces, where they evap- 
orate, to less curved ones through the bulk of the other 
phase. The corresponding growth exponent is z = 1/3 p|. 
At high temperatures (but still less than the critical tem- 
perature) the bulk diffusion is the only observed because 
it is faster than the other. Since bulk diffusion is a ther- 
mal activated process [Q, one expects that, lowering the 
temperature, a regime with z = 1/4 can be observed. The 
proposed form (||) for the mobility is able to catch these 
features. For shallow quenches (a <C 1) the mobility re- 
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mains constant in the whole system, while for very deep 
quenches (a = 1) r(<f) vanishes in the bulk phases where 
=1, suppressing the bulk diffusion. The diffusion along 
the interfaces is unaffected because r((j>) ~ 1 on domain 
boundaries. 

The effect of an order parameter-dependent mobility 
on systems with scalar order parameter has been stud- 
ied by simulations §,0]. It has been found that for a — 1 
the length scale L(t) grows as t 1 / 4 and for < a < 1 
there is a crossover between L(t) ~ i 1 / 4 and L(t) ~ t 1 / 3 . 
Recently, a non constant mobility has been also used to 
study the phase-ordering dynamics in systems with a vec- 
torial order parameter ||||. In H the expression (|[) for 
the mobility has been adopted. The limit N — > oo , where 
N is the number of vectorial components, is analytically 
solvable and corresponds to the self-consistent approxi- 
mation of the scalar case. It is characterized by a growth 
exponent 1/6 for the case a — 1. For < a < 1 the usual 
value 1 /4 of the growth exponent for vectorial systems is 
recovered in the asymptotic regime. A more general form 
of the mobility has been introduced in (^] , being given by 

r(0) = (i-0 2 r (3) 

where a is a positive real number. 

In this paper we are concerned with the phase sepa- 
ration of a binary mixture subject to an uniform shear 
flow when the mobility is given by Eq. (^) . When a shear 
flow is applied to a quenched binary mixture, the pat- 
tern of the phase-separating domains as well as the time 
evolution are strongly modified by the flow Q. Domains 
greatly elongated in the flow direction have been observed 



G. Gonnella el al: Phase-ordering dynamics of binary mixtures with field-dependent mobility in shear flow 



in simulations ]lC| , |ll[ and in experiments |12|] . A difference quantities is also calculated. Finally, in section 4 we in- 

Az = 0.8 -j- 1 between growth exponents has been mea- tegrate numerically the time evolution equation of the 

sured in some experiments, with the larger exponent in the structure factor and calculate the whole evolution of the 

flow direction [OJl4|. The deformations of domains affects physical observables as moments of the structure factor. 



the rheological properties of the system, giving rise to ex- In the asymptotic regime they are modulated by dumped 
cess stresses and to an increase Ar\ of the viscosity Jil|[l6|]. log-time oscillations. Our conclusions complete the article. 



2 The model 



This behavior has its explanation in the fact that mechan- 
ical energy is expended to deform the domains against 

the interfacial tension. When the domains are stretched 

n , , , , , , , . , , , , , We consider a binary mixture whose evolution is described 

to such an extent that they start to burst, the stored me- J 

, ■ i .... , , , . ., , by the diffusion-convection equation 

chanical energy is dissipated and the excess viscosity de- 

dH[4>(T,t)] 



creases. The phase separation in steady shear for binary d(f)(r, t) ^ ^ t)y) — V 
mixtures with constant mobility has been recently stud- 



r(0)v 



(4) 



ied in 17 , where the existence of an anisotropic dynami- 
cal scaling theory with different growth exponents in flow 
and in other directions has been shown. It is found that 
the excess viscosity, after reaching a maximum, relaxes to 
zero, exhibiting log-time periodic oscillations. Also other 



where the field <f>(r,t) describes the concentration differ- 
ence between the two components of the mixture, r(<p) is 
the mobility, which depends on the order parameter as in 
Eq. (H) . The order parameter (f> is convected by an external 
velocity field M . We choose a planar shear flow with 



physical observables are modulated by such oscillations, v = <yye x (5) 

which can be related to a cyclical mechanism of storing 



and dissipation of elastic energy. Here we want to see the 
effects of a non constant mobility on this scenario. 



where 7 is the shear rate, assumed constant, and e x is 
a unit vector in the flow direction. The free energy func- 
tional is chosen to be of the standard d 4 form: 



H\4>]= dv 



(G) 



The outline of the paper is as follows. In section 2 we 
present the model and derive the equation of time evolu- 
tion for the structure factor in a self-consistent approxi- where we assume that all parameters have been rescaled 



-^ 2 (r,t) + ^ 4 (r,t) + i|V0(r,t)| 2 



mation. In section 3 we report the asymptotic behavior of into dimensionless units [la ] and the system is below the 

the model. We found two different growth exponents for critical temperature. The two minima of the polynomial 

the flow and the shear directions, given in the leading seal- part of H[cj)] describe the pure states of the mixture, 
ing regime by z x = (5 + 2a)/2(2 + a) and z y = 1/2(2 + 0;), In this article we deal with the non-linear term of 

respectively. The asymptotic behavior of the rheological Eq. (0) in a self-consistent approximation; hence the term 
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(j) 3 appearing in the functional derivative 5H/6<j> is Lin- 3 The Scaling behavior 

earized as < <f> 2 > cf>, where < ... > stands for the average 

over the system. In the same way the term d> 2 in the mobil- » .., ,. „ ,, „ 

J Assuming simple scaling lor the structure factor, we write, 

ity is substituted by < <j> 2 >. In the Fourier space Eq. (§ for arbitrary gpace dimensionality dj 
becomes: 

d 

„ C(k,t) = JJJ2 i (t)F(X,r(7*)) (12) 

- (1 - S(t)) a k 2 [S(t) - 1 + fc 2 ] <£(k,t) (7) where the subscript i labels the space directions with i = 1 

along the flow, i?i is the average size of domains in the i-th 

where S(t) =< </> 2 (r,i) > and <p(k, i) is the Fourier trans- 
direction, X is a vector of components Xj = fci-Rj(t) and 

form of <A(r,t). 

F is a scaling function n. We also allow an explicit time 

The statistical quantity of experimental interest is the 

dependence of the structure factor through r(7t); notice 

time-dependent structure factor C (k, t) which is defined 

that since C(k, t) scales as the domains volume below the 

as < 0(k, t)<j)(— k, t) >. It obeys to the evolution equation: 

critical temperature, r must not introduce any further al- 
<9C(k, t) _ , 9 , . gebraic time dependence in C(k, t). From the numerical 

„ r „ n results of the next section, we will sec that F is a dumped 

- 2(1 - S(t)) a k 2 [S(t) - 1 + k 2 ] C(k, t) (8) 

periodic function of r. Inserting the form ( |l2| ) of C(k, i) 

which is closed by the self-consistency condition 



Our self-consistent approximation is equivalent to that 

— f ^ (Jfe f\ (9) made in the large- iV limit for vectorial systems pQ] . In this 

J\k\<q (27r) d i — i 

case simple scaling is not verified for 7 = when a — |21{] 

where q is a phenomenological cutoff. Rheological quanti- and a ^ g so that C ( kit ) has not the form (ffj) but can be 
ties of interest can be calculated as momentum integrals written as L d ^ {k/km) , where k m is the position of the maximum 
of the structure factor. Since we restrict the solution of the in the structure factor and £ is a function which depends con- 
model to the two-dimensional case, we consider the excess tinuously on k (multiscaling) . In the case with shear the exact 
viscosity A-q and the first normal stress ANi defined by solution is still lacking, so we cannot, in principle, rule out mul- 
[ P"q| tiscaling. However, since simple scaling is the leading approxi- 
1 f <ik mation in the regions of the maxima of the structure factor, we 

^=-7- l y k ^k X k y C(k,t) (10) 

' k ' <<? can use it to obtain the correct value of the growth exponents 

a , r f dk r 2 ,21 /i/i ,\ (apart from logarithmic corrections) because the momentum 

integrals which define the observables are dominated by the 
maxima of C(k,t). 
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into Eq. (0) we obtain 



dF 



1 X 1 F 2 = R 1 R^ 1 \t — 



i=i 



R-i 1 Ri(F + XiFi 



i - s(t)] \ 2 xf ( Rl 2 *l -i+s{t) ) F 



fc=l 



K13) 



where Fi = dF/dXi and a dot means a time derivative. 
Under the assumptions that R\ ^> Ri (i = 2, d) and 
Ri ~ R (i = 2, d) we can write 

~ f (5i^ 

7X1^2 = R~ 1 rA t — + R^R^F + X 1 F 1 ) 



r-^y: 



(F + X.F,) + 21- S(t) 



R-±xHj2x 2 k -(l-S(t))R 2 )F 



k=2 



(14) 



Since the l.h.s. of Eq.(|Ij) has no explicit algebraic time 
dependence, one has the asymptotic solutions 

Rl(t) ^ 7i (5+2«)/2(2+«) 
R(t) ~ tV2(2+«) 

(1 - 5(t)) ~ i 1 /^") (15) 
r( 7 £) - log 7* 

The growth exponents in the flow and in the shear direc- 
tions are z x = (5 + 2a)/2(2 + a) and z y = 1/2(2 + a). 

4 + 2a 

We observe that z x — 2„ = — j = 1. The value z v 

v 2(2 + a) y 

is the same found in Q, in the leading approximation, in 
a model having a field-dependent mobility with vectorial 
conserved order parameter without shear. Increasing the 
values of a, one obtains values smaller with respect to the 
case with constant mobility, when z x = 5/4 and z y — 1/4. 
The shear affects only the growth exponent z x which re- 
mains greater than 1 for every real and positive value of 



The previous arguments can be used to establish the 
scaling properties of the rheological coefficients. Inserting 
the form into Eq. ( fTo| ) we obtain 

Arj(t) ~ ( 7 t)-( 3 + a )/( 2 + Q ) 7 -( 1+a )/( 2 + Q ) 

x J XiX 2 F[X,T(t)]dX (16) 

Therefore, in the scaling regime, for each value of ft, the 
functions Ar\ corresponding to different values of 7 col- 
lapse each on the others if rescaled as Arj — > 7 ( 1+Q )/( 2 + Q )z\jy. 
A similar analysis can be done for the normal stress. It is 
straightforward to show that in the asymptotic regime 

AN! ~ ^!/(2+«) j X 2p [ X; dX 

Setting a — 0, we recover the previous results for the case 
with constant mobility [fl7| . 

4 Results and discussion 

In this section we consider the numerical solution of Eq. (|7|) 
We will present the results for the calculation of the aver- 
age size of domains and of the rheological indicators Aq 
and AN%, stressing the effects of a non-costant mobility. 
We solved the equation of time evolution of the struc- 
ture factor numerically in two dimensions, implementing 
a first-order Euler scheme with an adaptive mesh. The 
starting configuration for the structure factor is a con- 
stant value, which corresponds to a disordered state at 
very high temperature. Eq. (Q) has been solved for differ- 
ent values of 7. In the following, results will be shown for 
the case 7 = 0.01. Similar results have been obtained in 
the other cases. The two values a — 1 and a — 2 for the 



6 



G. Gonnella et al: Phase-ordering dynamics of binary mixtures with field-dependent mobility in shear flow 



Jt = 1 



jt = 4 



jt = 9 




\ J 




\ 




Fig. 1. The structure factor at consecutive times for 7 = 0.01 and a = 1. The k x coordinate assumes positive values on the 
right of the picture, while the k y is positive towards the upper part of the plane. The maximum value of \k y \ is 0.8. In the other 
direction, k x varies in the range [—0.8, 0.8], [—0.4, 0.4] and [—0.2, 0.2], respectively. The positions of the peaks on the left foil of 
the structure factor is (—0.22,0.37) for the highest and (—0.08, —0.24) for the other at ft — 4. At jt = 9 the highest peak is 
located at (-0.05, -0.23) and the other, at (-0.10, 0.36). 



parameter a appearing in @ have been considered. At the 
beginning the function C(k, t) develops a circular volcano 
shaped structure. This is then deformed, as consequence 
of shear, into an elliptic structure. The sizes of the axes of 
the ellipse decrease in time at a different rate, this being 
larger in the fc^-direction. During this evolution two dips 
start to develop in the volcano edge until C (k, t) is made 
of two foils. In Fig. 1 at jt = 1, the shape of C(k, t) rep- 
resentative of this stage of the evolution is plotted. Later, 
on each foil, two well-formed peaks can be seen. At jt ~ 4 
the peaks characterized by the larger values of \k y \ pre- 
vail. Observe that the structure factor is symmetric with 
respect to the change k — > k. The peaks with the smaller 
\k y I corresponding to a more isotropic configuration of do- 
mains, grow faster than the others until they prevail as it 
can be seen in Fig. 1 at yt = 9. These peaks continue 
to prevail along all the time evolution. What we observed 
is that the peaks continue always to grow and that the 



difference of their heights as a function of the shear strain 
jt is an increasing function being modulated by dumped 
oscillations. 

In order to have information about the growth of do- 
mains, we computed the typical domain size as 



Rx(t) 



J dkC(k, t) 
fdkk*C(k,t) 



1/2 



(18) 



and the same for the other direction. The values of R x and 
R y are plotted in Fig. 2 as function of the shear strain 
jt for a — 1 and a = 2. We plotted also the values for 
the case with constant mobility (a — 0). Some comments 
are in order here. The asymptotic behavior is the one ex- 
pected through the previous scaling analysis: for a = 1 
one has R x - t 7 ^ 6 and R y - t 1 / 6 ; for a = 2, R x ~ i 9 / 8 
and R y ~ i 1 / 8 . The growth exponents in both the direc- 
tions are decreasing functions of the mobility exponent a. 
This is reasonable since the mobility (||) becomes smaller 
and smaller when a increases, if (1 — <j) 2 ) < 1. Another 
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1 10 10 10' 

yt 

Fig. 2. The typical size of domains as function of the strain 
yt in the flow and in the shear directions for different values of 
a: a — (full line), a — 1 (dashed line), a = 2 (dotted line). 

consideration is about the superimposed log-time periodic 
oscillations. In the case with constant mobility these oscil- 
lations have an apparently constant amplitude. For non- 
zero values of a these oscillations are dumped. We will see 
that this feature is common to all the observables and will 
be discussed later. 

We turn now to the study of the rheological behavior 
of the system. The external velocity field causes additional 
stresses on the mixture. Important indicators are the ex- 
cess viscosity and the first normal stress. We calculated 
numerically Ar] and AN% using their definitions ( p^| ) and 
(pd|), respectively. The results are shown in Fig. 3 and 
Fig. 4. The excess viscosity reaches its maximum at the 
onset of the scaling when the domains are expected to be 
maximally stretched in the flow direction and the struc- 




1 10 10 10' 

yt 

Fig. 3. Plots of the excess viscosity vs the strain yt for different 
values of a: a — (full line), a = 1 (dashed line), a — 2 (dotted 
line). The inset shows the maxima of Arj as function of 7 for 
the case a = 1. The slope of the straight line is 0.8. 



10 - 



1 - 




10 1 o z 10 

Fig. 4. The first normal stress as function of the strain yt for 
different values of a: a — (full line), a = 1 (dashed line), 
a = 2 (dotted line). 
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ture of C(k, t) is the one shown at yt = 4 in Fig. 1. Ac- 
cording to our analysis of the scaling behavior, we expect 
the excess viscosity to scale with 7 as Ar\ ~ j-( 1 + a )/( 2 + a ) 
for fixed value of "ft. In the inset of Fig. 3 we report for 
the case a = 1 the dependance of At\m on 7. We find an 
exponent 0.8 slightly larger than the expected 2/3. The 
reason can be related to the fact that the excess viscos- 
ity reaches the maximum before the asymptotic regime is 
fully realized. Then Arj decreases as a consequence of the 
dissipation of the elastic energy stored by domains which 
start to burst when they are stretched furtherly. Therefore 
more isotropic patterns form and the typical structure of 
the function C(k, t) is the one at — 9 in Fig. 1. The ex- 
cess viscosity decreases with a power law behavior which is 
consistent with the predicted exponent — (3 + a)/(2 + a). 
The first normal stress reported in Fig. 4 decreases in 
time according to the exponent —1/(2 + a) after reaching 
a maximum. The amplitudes of all quantites plotted in 
Figures 2, 3, 4 are modulated by dumped log-time oscil- 
lations |22]j . We believe that the physical explanation for 
the dumping of oscillations may be found in the vanish- 
ing value of the mobility at equilibrium. In the case with 
constant mobility the origin of the oscillations is related 
to a cyclical mechanism of elongation and bursting of do- 
mains, which allows to store and dissipate elastic energy 
in the system [ ft7|| . In the present case, during the time 
evolution, the decreasing values of -T(</>) suppress diffu- 
sion in the bulk phase and inhibit the growth of small 
bubbles coming from bursting. Therefore, they cannot be 
stretched too much by the flow. In this way it is more dif- 
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ficult to store elastic energy in the system and the excess 
viscosity can increase only of a small amount. This mecha- 
nism of growth inihibition becomes stronger and stronger 
in the course of evolution causing the observed dumping 
of oscillations. 

In conclusion, we have studied the phase separation of 
a binary mixture with field-dependent mobility in shear 
flow. We proved that dynamical scaling holds for this sys- 
tem. There are different growth exponents in the flow and 
in the shear directions which depend on the mobility ex- 
ponent a. The difference in growth exponents is always 
1. All the physical observables have amplitudes decorated 
by dumped oscillations which are periodic in logarithmic 
time. We made a guess about the origin of this behavior. 
It would be an important endeavour to study this system 
in direct simulation of Eq. (Q) to deeply understand this 
phenomenon. 

We thank Federico Corberi for valuable discussions about the 
subject of this work. 
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